function h5imshow(h5file)

%%
h5 = hdf5info(h5file);

X = hdf5read(h5.GroupHierarchy.Datasets(1));
Y = hdf5read(h5.GroupHierarchy.Datasets(2));
% size(X)
% size(Y)
%%
% close all
imshows(squeeze(X(:,:,:,1:20)), 2);
imshows(squeeze(Y(:,:,:,1:20)), 2);
% imshows(cat(3, squeeze(X(:,:,:,1:5)), squeeze(Y(:,:,:,1:5))), 2);
% imshows(squeeze(Y(:,:,:,1:10)-X(:,:,:,1:10)), 2);
return

%%
% close all
% a = squeeze(X(:,:,:,1));
% b = squeeze(Y(:,:,:,1));
% figure; hist(a(:), 60); axis tight
% figure; hist(b(:), 60); axis tight
% figure; hist(b(:)-a(:), 60); axis tight